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Abstract: Ghost-free bimetric theories describe nonhnear interactions of massive and mass- 
less spin- 2 fields and, hence, provide a natural framework for investigating the phenomenon 
of partial masslessness for massive spin-2 fields at the nonlinear level. In this paper we ana- 
lyze the spectrum of the ghost-free bimetric theory in arbitrary dimensions. Using a recently 
proposed construction, we identify the candidate nonlinear partially massless (PM) theories. 
It is shown that, in a 2-derivative setup, nonlinear PM theories can exist only in 3 and 4 
dimensions. But on adding Lanczos-Lovelock terms to the bimetric action it is found that 
higher derivative nonlinear PM theories also exist in higher dimensions. This is consistent 
with existing results on the direct construction of cubic vertices with PM gauge symmetry. 
We obtain the candidate nonlinear PM theories in 5, 6 and 8 dimensions but show that none 
exist in 7 dimensions. 
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1 Introduction 

To focus attention and fix conventions, we start with a review of partial masslessness in linear 
Fierz-Pauli (FP) theory and the question of its nonlinear extension. We then describe our 
approach to the problem and summarize our results. 

1.1 Partial masslessness in Fierz-Pauli theory and beyond 

The linear fluctuations h^i, of a massive spin-2 field are governed by the Fierz-Pauli equation 
[1, 2] which in a d-dimensional spacetime with a background metric becomes, 

SjZ V - (V - ^G^^G^" V) + % (V - G^^G^" V) = . (1.1) 

A is the cosmological constant and rfipp is the mass. For mpp = 0, (1.1) must reduce to the 
linearized Einstein equation for a massless spin-2 field, hence £ is given by 



E^'^h — -- 



V- (1-2) 
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The structure of (1.1) is determined by the requirement of the absence of ghost. Generically, 
it allows for d{d — l)/2 — 1 propagating modes corresponding to a massive spin-2 field in 
d dimensions. The nature and dynamics of Gfj^u cannot be specified further in this framework, 
but a nonlinear setup that embeds (1.1) must address this issue. 

If G^u is a dS or, in general, an Einstein spacetime, then on the Higuchi bound [3, 4], 

m2p = ^A, (1.3) 

equation (1.1) develops a gauge invariance /i^j, — > h^^, + A/i^^ with [5-10], 

A V - (v.V. + e(x) . (1.4) 

This can be used to gauge away the helicity zero component of /i^^, so in d = 4 only the four 
polarizations ±2, ±1 survive. This is the linear "partially massless" (PM) theory.^ 

The obvious question of course is if the gauge symmetry associated with partial massless- 
ness can be generalized from the FP theory in dS spacetimes to a nonlinear theory of spin-2 
fields. Several recent studies have attempted to address this issue using different approaches 
[12, 14-16]. In particular, the authors in [12, 14] directly construct cubic vertices for h^^, with 
the above PM gauge invariance. This constructive approach makes interesting predictions 
about the nonlinear PM theory, showing that : 

• Cubic hfj^u interactions with a PM gauge invariance (1.4) exist only in 3 [12] and 4 
[12, 14] dimensions, as long as the theory involves no more than 2 derivatives. Hence, 
while linear PM theory exists in any dimension, 2-derivative nonlinear PM theories can 
exist only in 3 and 4 dimensions. 

• When higher derivative interactions are allowed, then PM gauge invariant cubic terms 
can be constructed even for dimensions d > 4 [12]. The structure of the higher derivative 
terms is such that for d = 4 one again recovers the 2-derivative theory. Hence the higher 
derivative terms are relevant only for d > 4. 

Any method of constructing nonlinear PM theories must give rise to the above features implied 
by the explicit cubic vertex calculations. 

In [17] we obtained a potential nonlinear PM theory as a special ghost-free bimetric 
theory [18] in d = 4. There it was not obvious that this construction was consistent with 
the above dimension dependent features of PM theories. In the present paper, we implement 
this construction in arbitrary dimensions obtaining candidate PM theories including higher 
derivative terms. It is found that the construction indeed meets the expectations from the 
cubic vertex analysis. Let us consider this in some detail. 

It is natural to look for nonlinear PM theories among the nonlinear extensions of the 
FP theory (1.1) that do not suffer from the Boulware-Deser ghost instability [19, 20]. For 

^PM theories also arise in the more general context of higher-spin theories, see for example, [11-13]. But 
here we concentrate only on the spin-2 case. 
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Gfj,u = Vfiu such a nonlinear theory was obtained and shown to be ghost-free in a certain 
"decouphng hmit" in [21, 22] (also see [23]), whereas the complete nonlinear proof of absence 
of ghost was given in [24, 25]. However, these models do not admit dS solutions with massive 
FP perturbations around them. A generalization with extra parameters that accommodates 
a non-flat but non-dynamical Gf_iu was first considered in [26]. In this case decoupling limit 
arguments do not exist to argue the absence of the BD ghost, but unitarity was proven 
directly in the nonlinear theory in a Hamiltonian analysis [25, 27, 28]. Finally, the completely 
dynamical ghost-free theory was given in [18], as a bimetric theory conveniently formulated 
in terms of two interacting spin-2 fields Qfj^i, and /^jjy. This theory contains the two constraints 
required to eliminate the BD ghost [18, 25]. For related work see [29-48] 

As such, neither of the two metrics in the bimetric theory corresponds directly to the 
massive fluctuation in (1.1). But the theory can be reformulated in terms of a massive spin-2 
field coupled to a massless metric G^i, such that linear perturbations (^M^ around a 
background = always satisfy the FP equation (1.1) where now one identifies /i^^ = 
(5M^ [49]. Furthermore, (1.1) is also supplemented by a massless equation for SG^y, while 
G^y is determined as a solution to an ordinary Einstein's equation. In the absence of matter 
couplings, a background = forces G^y to be an Einstein metric.^ Hence, the bimetric 
theory provides a completely dynamical and ghost-free nonlinear extension of (1.1). For this 
reason it also provides the natural arena for investigating partial masslessness beyond the FP 
theory. 

Bimetric theories contain several parameters that can be easily tuned to put equation 
(1.1) in the dS backgrounds on the Higuchi bound (1.3). The linear theory will now have the 
gauge symmetry (1.4), but only the ^ = constant part of this transformation preserves the dS 
backgrounds. In [17] it was argued that the consistency of these dS preserving transformations 
with the dynamical nature of the backgrounds leads to a criterion for partial masslessness that 
is powerful enough to fix most of the remaining parameters of the model. The criterion was 
implemented m. d = A leading to a unique class of nonlinear candidate PM theories [17]. The 
complete gauge symmetry of the nonlinear model is not yet known but some special cases can 
be considered. Here we implement this criterion in arbitrary dimensions to obtain candidate 
PM theories that are consistent with expectations from the cubic vertex calculations described 
above. 

1.2 Summary of results 

Before getting into technical details, let us briefly summarize our results. 

• The bimetric theory has a well defined mass spectrum around proportional backgrounds 
fpLv = c^g^Lv where c is generically determined by the equations of motion. We obtain 
the mass spectrum of the fluctuations in arbitrary dimensions with explicit expressions 
for the mass and the cosmological constant. We also construct the nonlinear extensions 

^The theory has a much more compUcated behaviour around backgrounds with M^^ 7^ but this is not 
relevant for comparison to the FP equation. 
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of the massless and massive fluctuations which are useful for the purpose of fixing the 
normalizations in the PM case. 

• We then discuss the linear PM theory in the context of bimetric theory linearized around 
dS backgrounds. The consistency of a subset of (1.4) with the dynamical backgrounds 
leads to a simple criterion for identifying nonlinear PM theories as special bimetric 
theories with the parameters fixed such that the equations of motion leave the c in 
the ansatz f^^i, = c^g^u undetermined. In the PM theory c can be traded off with a 
constant gauge parameter. Applying this criterion to the 2-derivative bimetric action 
we find that nonlinear PM theories can exist only in c? = 3 and d = 4. Although the 
full gauge transformation of nonlinear PM theories is not yet known, we show that 
the mass, cosmological constant, couplings and the background metric Gfj^u all become 
c-independent, and hence gauge invariant, for the PM parameter values. 

• To explore the possibility of recovering nonlinear PM theories for d> A with the help of 
higher derivative terms, we consider the bimetric theory with extra Lanczos-Lovelock 
terms for both metrics.^ The analysis of the spectrum around dS backgrounds can easily 
generalized to the LL terms using recent results of [51] obtained for the LL extension 
of the Einstein-Hilbert gravity. We show that the construction outlined above can now 
be carried out leading to potential nonlinear PM theories even for d > 4. In particular, 
we obtain such nonlinear theory for d = 5, 6, 8 but show that no PM theory exists for 
d = 7. 

• We also outline the structure of the bimetric action expressed in terms of the nonlinear 
extensions of the massless and the massive fields. This will be useful for comparison 
with the outcome of the direct cubic vertex construction, although such an explicit 
comparison has not been attempted here. 

2 Bimetric description of massless and massive spin-2 fields in d dimensions 

In this section we consider the ghost-free bimetric theory in d dimensions and describe in what 
sense it is a theory of a massive and a massless spin-2 fields. To this end, we consider a class 
of bimetric backgrounds, the proportional backgrounds, around which linear perturbations 
decompose into well defined massless and massive modes. These can then be promoted to 
nonlinear fields, generalizing the notion of massless and massive spin-2 fields beyond the 
special class of proportional backgrounds. The analysis generalizes the d = A case considered 
in [49]. 



Such an extension has been considered earlier in a different context in [50]. 
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2.1 The bimetric theory in d dimensions 

The ghost-free bimetric action for two spin-2 field Qfj^u and /^^ [18] can be easily generalized 
to arbitrary dimensions,^ 

Sgf = j d^x [m1-^^MR{9) + mf^^ATW) " 2m''^\V{S-l3n)\ , (2.1) 

where \/\g\ = \/\ det g\ and S simply stands for the matrix square-root, 

S^^/^. (2.2) 

The interaction potential is given by, 

d 

V{S;(3n) = Y,Pnen{S), (2.3) 

n=0 

where en{S) are the elementary symmetric polynomials of eigenvalues of S. They are ex- 
pressible as polynomials of Tr(S'''') and can be iteratively constructed starting with 60(5*) = 1 
and using Newton's identities, 

n 

en{S) = — y^i-l f Tt{S'') en-k{S) . (2.4) 

k=l 

In particular, 6^(5) = det S and e„(S') = for n > d. 

The in (2.3) and the Planck masses nig and mf are the d + 3 free parameters of the 
theory. The mass parameter m is degenerate with the overall scale of the /?„. The action 
contains the cosmological terms /3o\/M" and I3d\/\f\, but the actual cosmological constants 
are to be read off from the equations of motion. 

The absence of the Boulware-Deser ghost in bimetric theory in d = 4 was proved in 
[18, 25]. It is straightforward to extend the ghost analysis to arbitrary dimensions, especially 
using the "deformed determinant" representation for the potential [26] . Setting to a non- 
dynamical flat metric and tuning the to exclude a cosmological constant contribution in a 
flat Qfj^i, background, one recovers the massive gravity model of [22]. 

The bimetric potential (2.3) has the following useful symmetry property under the inter- 
change of 5^iy and fij,i, that sends S ^ [18], 

y^\ViS■,|3n) = ^AT\V{S-'■,Pd-n). (2.5) 

This allows us to obtain the /-sector equations from the g-sector of the theory. 

The sourceless equations of motion obtained on varying the action (2.1) with respect to 
g^"' and ff"" are, 

R,M - b>..R{9) + = , R^M) - y,uR{f) + ;^V^, = 0. (2.6) 



^In the paper we use the sign and curvature conventions of Wald [52] so that i?^^ — daTf^i, + ■ 
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The potential contributions have the form, 

= 9,x Vl'^iS) , y/, = Vl^iS-') , (2.7) 

and are exphcitly given by (A. 8) and (A. 9) in the the appendix (where the details of the 
derivation are also provided). 

2.2 The proportional backgrounds 

Let us now consider a class of solutions to the bimetric equations in which the metrics are 
proportional to each other, 

= c'5mv (2.8) 

Such backgrounds have two important properties: 1) they are the most general class of back- 
grounds around which the bimetric theory has well defined massless and massive fluctuations 
with a Fierz-Pauli mass term, 2) they coincide with classical solutions in general relativity.^ 
Indeed, for this ansatz, the bimetric equations (2.6) imply that c is constant and then reduce 
to two copies of cosmological Einstein's equation for (^^j/, 

R^L^g) - \9i.uR{g) + 5MvAg = , R^,^{g) - \gf,yR{g) + g^u^f = . (2.9) 
The cosmological constants are given by (for details see the appendix), 

y n=0 / n=l 

where, (^) = is the combinatorial factor. Note that does not contain fid while A^- 

is independent of /3o- 

The consistency of the two equations in (2.9) with each other then implies,^ 

Kg=Kf. (2.11) 

This is a polynomial equation that, generically, determines c in terms of the parameters of 
the theory, the exception being the PM case [17] to be discussed later. In particular, the 
background equations (2.9) admit de Sitter solutions that are relevant for the identification 
of the linear PM theory. 

Since (2.11) is homogeneous in the /3fc, the c determined by it is independent of the overall 
scale of the /3fc and will depend on at most d + \ parameters, say /3i//3o, • • • Pd/ f^o and a = 
nif/mg. We are interested in parameter values that result in real non-zero c, as c also appears 
in other quantities, like the Fierz-Pauli mass of the fluctuations. In practice, such parameter 
ranges can be easily found since, given any d of the parameters, say, {a,/32/f3o, ■ ■ ■ ,f3d/Po}, 
one can determine the allowed range of the remaining one, /3i//3o, by expressing it as a function 
of the real c using (2.11). 

^For another feature of such solutions, see [53] 

®The discussion can be easily extended to include sources, in which case, (2.8) also forces a proportionality 
relation between the two energy-momentum tensors as discussed in [49] . 
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2.3 Linear mass eigenstates 

To determine the mass spectrum of the theory, let us now consider hnear perturbations around 
the proportional backgrounds, 

QflU = (jflU + Sg^u , ffXV = Uv + 5f^y , (2.12) 

with f^u = c^g^u- It turns out that only around such backgrounds the fluctuations can be 
combined into definite mass eigenstates. A massive mode in a general background is identified 
through the appearance of a Fierz-Pauli mass term. On linearizing the equations of motion 
(2.6), one gets (for details, see the appendix) 

^'^u^dpa - 3^A, (^6g^, - ^g^ug^'Sgp^^ + ;^Ng^x (tv{SS)6^ - 6Si) = , (2.13) 

^P^fp^ - dh^9 ('^^^ - ^5m.5''"5/p<x) - ^Nc^-^g^x (Tr{6S)6t - 5Si) =0. (2.14) 

S^u is defined through (1.2), now with background metric g^u, and N is given by (A. 26). 

= ^g^''^{Sf — (?^g)\u enters both equations in the FP combination, hence we expect to 
get a massive mode JM^jy ~ g^i\^S\. 

The linearized equations are easily diagonalized in terms of a massless fluctuation JG^^ 



and a massive fluctuation 6Mfj_u, 



6G^, = A{c) (^%, + c^-*a''-'<5/^.j , (2.15) 
6M^, = ^ {6f^, - c^Sg^,) . (2.16) 
A{c) and B{c) are normalizations to be determined later and we have used the notation, 

a = ^. (2.17) 

nig 

Indeed, adding (2.13) and (2.14) to cancel the FP mass term gives a massless spin-2 equation, 

£pSGp^ - ^Ag (^6G^, - ^g^.r^SGp^^ = . (2.18) 

On the other hand, subtracting the right combination of (2.13) and (2.14) gives the FP 
equation (1.1) for the massive spin-2 fluctuation, 

2 m? 
6Mp^ - j^^g {5M^, - Ig^.gP'^dMp^) + {5M^, - g^,,gP"8Mp^) = . (2.19) 

Here, the FP mass is given by, 

2 m<^ /l + (ac)'^-2x ^ /^_2\ 



We emphasize that equations (2.18) (2.19) are written with g^y as the background metric. 
Then the expressions for mpp (2.20) and Kg (2.10) refer to this background metric choice. 
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2.4 The nonlinear massive and massless fields 



The massless and massive fluctuations can be regarded as perturbations of some nonlinear 
fields G^jy and Af^. If this is the only criterion, the choice of nonlinear fields is far from 
unique. However, we also require that the relation between the nonlinear fields and the 
original variables g^i^ and /^^ is simple enough that it is invertible in a useful way [49] . These 
criteria single out the following straightforward nonlinear extensions of the mass eigenstates, 



G^, = A{c) [g^, + c'^-^a'^^'/^.j , (2.21) 
M^, = B{c) {G^pSP, - cG^,) . (2.22) 



where S = \l g~^f and the dimension dependent normalizations A and B will be fixed later 
in the context of the PM theory. It is not claimed that G and M*-^ propagate respectively 2 
and 5 degrees of freedom nonlinearly. 

The bimetric action (2.1) can be re-expressed in terms of the new nonlinear fields as a 
theory of a massive spin-2 field interacting with a massless spin-2 field G^,y. However, it 
turns out that Gfj_,y cannot be coupled to matter in the standard way without reintroducing 
the Boulware-Deser ghost and, hence, it cannot be regarded as the physical gravitational 
metric [49]. In other words, the spin-2 mass eigenstates differ from the states produced by 
bimetric interactions, as is familiar from other contexts in particle physics.^ However, the 
theory in terms of G and M'^ is relevant for discussing partial masslessness which is known 
to arise only in the absence of matter couplings. 

On the proportional backgrounds, the nonlinear fields reduce to, 

d-2\ - r^G 



Gp, = A{c) (^1 + {acf-') g^u , = . (2.23) 

A non-vanishing expectation value for signals non-proportional backgrounds and hence 
parametrizes deviations of the bimetric solutions from solutions in general relativity. 

The linearized equations (2.18), (2.19) are written with g^i, as the background metric. 
To describe the theory entirely in terms of the massless and massive modes, we rewrite these 
in terms of G^jy using (2.23), 

SP SGp^ - {^G^, - \Gp,&''bGp,) = , (2.24) 

SP^, 5Mp„ - j^Ag {5M^, - ^G^,GP^5Mp^) + {5M^, - Gp,&^5Mp„) = . (2.25) 

Here Sjiu is given by (1.2). In the new background convention, the mass and the cosmological 
constant are rescaled to, 

" A{c) (1 + (ac)'^-2) ' '^9 - ^(c) (1 + (ac)'i"2) • l^-^^^ 



^To describe a massive spin-2 field interacting with a gravitational metric that is also sourced by conventional 
matter, one should use the fields gip^ and M^^ = g^AS^ — cgfj,„ [49]. 
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3 Partial masslessness and Bimetric theory in arbitrary dimensions 

To summarize, on linearizing, the bimetric action leads to the FP equation for massive spin-2 
fields in a cosmological background (2.25), along with equation (2.24) for massless fluctuations 
of the background itself. Hence it provides a completely dynamical nonlinear generalization 
of the FP equation discussed in section 1.1. In particular, in this setup, one can investigate 
partial masslessness beyond the linear level. 

Working in 4-dimensions, ref. [17] obtained a criterion to identify the nonlinear PM 
theory among the set of ghost-free bimetric theories. Here the considerations are extended to 
d-dimensions showing that, while linear PM theories exist for any d, nonlinear PM theories 
based on the bimetric action (2.1) exist only in d = 3 and d = A. This is consistent with the 
results of [12] which shows that, in a 2-derivative theory, spin-2 cubic interactions invariant 
under the linear PM gauge symmetry exist only in these dimensions. Furthermore, [12] 
shows that when the 2-derivative restriction is relaxed, cubic terms with PM symmetry can 
be constructed in higher dimensions as well. This case is considered in the next section. 

It should be emphasized that at this stage we do not identify the PM gauge symmetry in 
the nonlinear theory, except for a subset of the transformations around maximally symmetric 
backgrounds. Hence the claim is that if a nonlinear theory with a PM gauge symmetry exits, 
it must belong to the set of theories identified here. 

3.1 Linear PM symmetry in bimetric variables 

Since equation (2.25) coincides with the FP equation (1.1), it is straightforward to identify 
the linear PM theory and the associated gauge symmetry (1.4) in bimetric theory. Hence, 
when the Higuchi bound is satisfied, 

2 

mlp = j—^kg, (3.1) 

equation (2.25) develops a new gauge invariance that transforms the massive spin-2 fluctuation 
as (5M^jy — )• 5M^v + A{6Mfj_^), where 

A(5M^,) = (V^V, + ^^_|^k__G^^)e(^) . (3.2) 

Now we also have an equation (2.24) for the massless fluctuations SG^^- Since the massless 
equation has no unusual symmetries, the above transformation must be supplemented by, 

A{5G^,) = (3.3) 

Using (2.15) and (2.16), one obtains the corresponding transformations of 6g and 6f as. 

In any dimension, these equations realize the linear PM gauge transformations in the bimetric 
theory around proportional backgrounds. The associated Higuchi bound (3.1) simply reduces 
the number of bimetric parameters by just one and it too can be satisfied in any dimension. 
But we will see that the nonlinear extension is much more restrictive. 
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3.2 Identifying the nonlinear PM theory in d dimensions 

The gauge transformations (3.4) of the hnear PM theory lead to the criteria for identifying 
potential nonlinear PM theories provided one works in a setup that treats the backgrounds 
dynamically. The argument is as follows [17]: A nonlinear PM theory must be invariant under 
an extension of (3.4) that transforms the nonlinear fields as, say, 

9fj.u ~^ g^u ~ 9fj.u + ^g^jLu 1 ffiu ~^ ffj,u — f^J,v + ^ffiu ■ (3-5) 

For g = g + 6g, the transformed field also splits as g' = g' + 6g' , but there is a gauge ambiguity 
depending on how Ag is split. For instance, a small Ag may be viewed as transforming either 
the background, g' = g + Ag, or the fluctuation, 6g' = 5g + Ag. The same holds for /. 

The implication of this for the linearized bimetric theory is that the variation A{6g) of Sg 
(3.4) can, alternatively, be regarded as a transformation of the background g to g' = g+A{Sg), 
keeping 6g unchanged. The same holds for f^^. Let's adopt this latter point of view. Now, 
for a generic gauge parameter ,^(x) in (3.2), the new (g' , /') are not proportional backgrounds 
and hence are not dS metrics, in which case not much is known about partial masslessness. 
To keep to dS backgrounds we restrict A{5g) and A{6f) to constant ^(x) = and rename 
them to Ao(g), Ao(/) to emphasize that they transform {g,f) rather than {6g,6f). Then, 
given f = g, the new backgrounds are related through, 

f' = c\^o)g, (3.6) 
and the dS preserving gauge transformations at background level are, 

g' = g + Ao{g), /' = / + Ao(/), c' = c + Ao(c). (3.7) 
The Ao-variations above are obtained from (3.4), (3.2) and (3.6) for constant ^(x) = 

Ao{g) = --{acf~^gXd^o, Aoif) = 2cg X^Co , Ao(c) = fl + (qc)'^"^! ^^^0 , (3.8) 
c L J 

where, stands for, 

A 2Ag 

A and B are the normalizations of the nonlinear massless and massive fields Gfj,,y and to 
be determined later. 

Thus, the conclusion is that if a nonlinear theory with PM gauge symmetry exits, its 
equations of motion on proportional backgrounds must be invariant under the transformations 
(3.7). This requirement puts constraints on the bimetric theory (2.1) and singles out the 
potential nonlinear PM cases. The criteria for the invariance are listed below. 

1. Recall that for proportional backgrounds / = c^g, the bimetric equations of motion 
imply Ag = Af (2.11), explicitly, 

n=0 ^ ^ n=l ^ ^ 
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Generically, this fixes c = c{a, Pi) and excludes the possibility of invariance under trans- 
formations (3.7) that require changing c. Thus the nonlinear theory can be invariant 
under (3.7) only for special values of (/3„,a) for which (3.10) leaves c undetermined. 
This is the necessary condition for the existence of nonlinear PM theories for these 
(/3n, a) values. 



2. When the theory is expressed in term of the massive and massless fields and 



the FP mass and the cosmological constant take the form (2.26), or explicitly, 

'd-2' 



rhpp 



m 



d-2 



A{c) 



n=l 



n 



1 



A. 



m 



m„ 



d-2 



n=0 



d-1 



n 



(3.11) 



(3.12) 



Since the gauge transformations in the nonlinear PM theory involve shifts of c, then 
TTipp and Ag, as well as the effective couplings, must become c-independent. 

3. The gauge symmetry of the linear PM theory leaves the massless fluctuation 5G^i, 
invariant. Hence at the background level too, transformations (3.7) must keep the 
corresponding nonlinear background G^i, invariant. This requirement determines the 
normalization A{c) of G^^. Indeed from (2.23) and (3.7) one can compute. 



g/iuA 



1 + (ac 



4-2 



Ao(c)-^ln(^ 

ac 



1 + {acf 



0. 



which fixes the normalization, upto a c-independent factor, as. 



A{c) = 1 + (ac) 



d-2 



Then, the background massless and massive fields become, 



G 



1 + (ac) 



d-2 



M. 



0. 



(3.13) 



(3.14) 



(3.15) 



As discussed below, the first criterion is satisfiable only in 3 and 4 dimensions (not counting 
d < 3 which has trivial dynamics). Then in these cases, the second criterion is satisfied only 
for the normalizations A{c) determined by the third criterion. Hence, the potential nonlinear 
PM theories in 3 and 4 dimensions are manifestly invariant under the dS preserving part of 
the PM gauge transformations. Let us now consider these cases separately. 

3.3 Nonlinear PM theory for d=3, 4 

In 3 dimensions the condition (3.10) that generically determines the c in / = c^g reads 



(aPo - pi) + 2(a/3i - ^2)0 + (a/32 - Pz)c^ 



0. 



(3.16) 
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It leaves c undetermined only when the coefficients of all powers of c vanish, 



/3i = a/3o, /32 = af3i, = , (3.17) 

or simply, /3„ = a"/3o for n = 1,2,3. Only for these parameter values the restricted gauge 
transformations (3.7) lead to new proportional backgrounds /' = d'^g' with c' 7^ c, without 
violating the bimetric equations of motion. For the above parameter values, the mass and 
cosmological constant (3.11, (3.12) measured in the metric G^^ become. 



m 



3 



mfp = Ag = /3i . (3.18) 

nif 

Thus they satisfy the Higuchi bound and are indeed independent of c and hence gauge in- 
variant, but only for the normalization ^ = [1 + (ac)] given by (3.14). 

In 4 dimensions (first considered in [17]), the condition (3.10) becomes, 

/?i + (3/^2 - a^/3o) c + (3/33 - ^a'^^i) + {p4 - + a^li^c^ = . (3.19) 

This leaves c undetermined for 

/?i = /33 = , a^/3o = 3a2/32 = 15^ . (3.20) 

Then the nonlinear PM theory must correspond to this choice of parameters. For these 
parameters, the mass and cosmological constant satisfy the Higuchi bound, 

2 

-rnlp = -Ag = 2 — /32 . (3.21) 
3 rrij 

Again, they are independent of c, and hence gauge invariant, precisely for the normalization 
A = \ (3.14) that also renders G^u invariant. 

The bimetric action (2.1) with the /3's given above satisfies the necessary conditions for 
being a nonlinear PM theory for (i = 3, 4. The PM gauge symmetry of the nonlinear action 
(2.1) is not yet known. But the discussion shows that, at least for proportional backgrounds, 
the nonlinear equations are invariant under the dS preserving subset of the PM gauge trans- 
formations (3.7). In fact, this subset can be further extended to the nonlinear dS preserving 
transformations , 



d = c + a, g' 



2 

d-2 



^ + (3.22) 



1 + (Q(c-ha))'^-2 

where a is a finite parameter and the transformation of / follows from / = c^g. To linear 
order in a = Ao(c), this reduces to (3.7). The above transformations keep the background 
G^i, (2.23) unchanged and also leave the background g^i, and f^i, equations of motion (2.9) 
invariant. 
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3.4 The quadratic action 

As a last check, let us consider the quadratic action for the fluctuations 6Gf^u and dM^j^i, given 
by (2.15), (2.16) and (3.14) with B{c) = 1, that is. 



6G 



6M, 



1 + {acf 

Yc ^^^^^ " 



5gnv + C' 



(3.23) 
(3.24) 



This action is needed to determine the eff'ective Planck masses of 5Gfj,i, and 5M^y and to verify 
that they are indeed independent of c and hence gauge invariant. Expanding the bimetric 
action (2.1) to second order and expressing the result in terms of the above combinations one 
obtains, 



5(2) 



m„ 



,d-2 



6MS6M 



IG 



-5G85G + 



(d-2) 



Tr[(5M^ 



{d-2 
1 



Tr[(5G2 



\TT[5Gf 



Tt[5M\ 



+ 



m 



FP 



Tr[5i\/]2 - Tr[5M2 



(3.25) 



Here, rhpp and Ag are given by (3.11), (3.12) which become independent of c for the PM 
parameter values. Also note that the Planck masses are independent of c for the correct nor- 
malization in (3.23) determined by other considerations. Furthermore, the quadratic action 
fixes B(c) = 1 since otherwise, the 5M sector will have a c-dependent coupling. 

In the appendix we outline the structure of the complete nonlinear action expressed in 
terms of G^jy and M^. This is useful for computing the cubic and higher order interactions 
of the PM field to be compared with the explicit cubic vertex computations in [12, 14, 16]. 
We do not attempt such a comparison here. 

3.5 Absence of nonlinear PM theory in d>4 

The necessary condition for the existence of nonlinear PM theories is that (3.10) leaves c 
undetermined. On relabeling the sum on the right hand side, this equation can be recast as, 
/ d-l 2 \ 



a 



d-2 



n+d—2 I C 



0. 



(3.26) 



n=0 



71=3— d 



where, 



A-n 



(3.27) 



'd-V 
n 

To leave c undetermined, the coefficients of each power of c must vanish separately. For d > 4, 
this gives the following 3 sets of equations: 



3- d<n< -1 
< n < 2 
3<<n<d-l 



Pn+d~2 

a 

/3n = 0. 



0. 



d-2 A n 
-^n Hn 



B. 



n+d-2 



Pn+d- 



0, 



(3.28) 
(3.29) 
(3.30) 
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The first set implies that Pm = for 1 < m < d — 3 whereas the last set implies that /3m = 
for 3 < m < d — 1. Then the second set implies Po = Pd = 0. Hence, /3„ = for n = 0, • • • , d. 

To recapitulate, the bimetric action (2.1), when linearized around proportional back- 
grounds, leads to the massive FP equation (2.25). On the Higuchi bound (3.1), this gives a 
linear PM theory in any dimension. The consistency of the dS-preserving part of the linear 
PM gauge symmetry with dynamical backgrounds provided the necessary (bot not sufficient) 
condition for the existence of nonlinear PM theories. This criterion was satisfied in 3 and 4 
dimensions but not for d > 4. 



4 Partial Masslessness with Lanczos-Lovelock terms 

Our nonlinear results are consistent with the observation in [12] that in a 2-derivative theory, 
cubic spin-2 interactions with PM gauge symmetry exist only in 3 and 4 dimensions. However, 
[12] also finds that in the presence of higher derivative terms, cubic PM interactions exist even 
for d > 4. These higher derivative terms reduce to 2-derivative terms in d = 4. This hints that 
the nonlinear higher derivative interactions correspond to Lanczos-Lovelock (LL) terms which 
will also preserve the unitarity of the bimetric theory. This section considers the LL extension 
of the bimetric theory from the point of view of the PM analysis. Such a theory has been 
considered in [50] in a different context. Recently, in [51] the spectrum of Einstein-Hilbert 
gravity with LL terms was studied in dS backgrounds. Here, this analysis is generalized to 
the bimetric theory and used to investigate the existence of PM theories beyond d = 4. 

To summarize, it is shown that, for / = c^g and around dS backgrounds, the structure 
of massless and massive fiuctuations in bimetric theory with LL terms is exactly the same as 
in pure bimetric theory, though with modified parameters. Hence the PM analysis previously 
developed for the pure bimetric case also applies here. In particular, the linearized theory 
has PM gauge symmetry on the Higuchi bound. Then for a nonlinear extension of the theory 
to exist, the necessary condition is that the equations of motion must leave c undetermined. 
This condition is easily implemented if one obtains a polynomial equation for c. 

For d < 4 the LL terms are of no consequence and the pure bimetric analysis goes through. 
The first new contribution is from the quadratic LL term for d > 5 and this remains the only 
contribution for d = 5,6. In this case, one can easily obtain a polynomial equation for c and 
determine the potential PM theory. The cubic LL term starts contributing for d > 6 but 
in this case one does not automatically obtain a polynomial equation for c. However, using 
a decomposition into the Grobner basis one again obtains a polynomial equation for c, as 
discussed in subsection 4.3. The analytical derivation of the PM theory in d = 5 is presented 
in subsection 4.4, whereas the PM theories for d = 6, 7, 8 are determined with the help of a 
computer and are discussed in subsection 4.5. Beyond this, the analysis gets more involved 
and is not attempted here. 
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4.1 Bimetric theory with Lanczos-Lovelock terms in dS spacetimes 

Adding higher order LL invariants for both g^^, and /^^ to the bimetric action gives, 



n=2 

'x^v(^./^f;Pn) , (4.1 



[d/2] 
n=2 



where In and In are couphng constants of mass dimension — 2(n — 1). The sum terminates 
at [d/2] (integer part of) since the Lovelock terms vanish identically when d < 2n and are 
topological invariants when d = 2n. Our conventions for the Lovelock invariants are, 

n 

~ 2" "i/^i---""-^" J-J- z^'-'''- ' "i/3i---"n/3n ~ n! PI ' ' ' "" Pi] ■ ^ ' 
Note that the action (4.1) is invariant under simultaneous interchange of, 

^a'^V, /3„^a"/3d_„, ll^c?"''^ll. (4.3) 
The equations of motion obtained from this action read, 

[d/2] 

Ry^M - \9^vR{9) + ^n^£^(5) + -P^^A = 0' (4-4) 

n=2 
[d/2] 

i?^.(/) - ^9/..i?(/) + E 'n^i"^(/) + J^^/v = ' (4-5) 

n=2 ^ 

(n) 

where the Lovelock tensors Q^iu appear as the result of varying the Lovelock invariants. 

We are interested in the maximally symmetric solutions of the above equations that allow 
for massive spin-2 excitations. In the bimetric setup, these belong to the class of proportional 
backgrounds ff^i, = (?gf_iu- For a maximally symmetric spacetime with cosmological constant 
A, the curvatures are, 

R^iupa{g) = _ _ 2^ idtipQua - gup9i,a) , R^iuio) = JZ^diiu , Rio) = , (4.6) 

with the corresponding equations for f^y. With this ansatz it is enough to consider the traces 
of the equations of motion (4.4) and (4.5) to obtain two equations that generically determine 
A and c in terms of the parameters of the theory. The computation is simplified by, 

5^^ai:)(5) = ^^An)(ff), (4-7) 

which allows us to work with the Lovelock invariants, rather than the corresponding Lovelock 
tensors. Furthermore, for the proportional ansatz we have, 

>C(„)(/) = c-2-£(„)(5). (4.8) 
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Now, for the ansatz (4.6), the Lovelock invariants become. 



2" J' 

C,. = NJd)X^ , with NJd) = — — — . (4.9) 

^''> ^ ' ^ (d- l)"(d-2)"((i-2n)! ^ ^ 

Using the above relations, the traced equations of motion (4.4) and (4.5) give, 

f'/'' d-2n 

^+J2^n -^^n{dW - A. = , (4.10) 

n=2 

[dm , _ „ 

A + E ^'"'"^^ - A/ = . (4.11) 

n=2 

The contributions A^j from the bimetric potential are still given by (2.10). In general, these 
equations determine the cosmological constant A and the proportionality constant c. Hence, 
they specify the maximally symmetric solutions in bimetric theory with LL terms. 

Now, we consider the spectrum of fluctuations in the above dS backgrounds. The analysis 
is greatly simplified using results of [51], which considered Einstein-Hilbert gravity with LL 
terms, Seh+ll- For a single metric gfj,^ = gij.u + ^9^iu described by Seh+ll, [51] showed 
that the quadratic action for 5gfj^iy in an (A)dS background gfj^iy was exactly the same as the 
quadratic action in pure Einstein-Hilbert gravity in the same background but with a modified, 
effective, Planck mass. Applied to the bimetric action with LL terms (4.1) this implies that 
at the quadratic level, the theory has the same structure as the pure bimetric theory (2.1) 
discussed earlier, but now with modified Planck masses. 



- d~2 ^d~2 

iJij = rrij 



-(-3,iEf^(<,3^) c], 



n=2 



This implies that, just as before, the spectrum of fluctuations consists of massive and massless 
spin-2 modes in a dS background with cosmological constant A. In particular, the Fierz-Pauli 
mass has the same form as in (2.20), but now involves the modified Planck masses, 

m I / nif \ \ \ - I a — I 



k=l 



Since the massive mode satisfies the FP equation, it is obvious that the linear theory will 
exhibit partial masslessness at the Higuchi bound, m|p = ^jrrA- Below we show that in the 
presence of the LL terms, nonlinear potentially PM theories can be found even for d > 4. 
Before this, let us look at the quadratic and quartic LL terms in more detail. 



- 16 - 



4.2 The quadratic Lanczos-Lovelock term 

We now apply the general results of the previous section to the quadratic LL term which is 
the Gauss-Bonnet combination, 

£(2) = R^- AR^^uR^" + Rf^upaR^""'' . (4.15) 

This contribution is trivial for d < 4 but modifies the bimetric theory for d > 4. For a 
maximally symmetric spacetime with cosmological constant A it reduces to 



{d-l){d-2)' 



^2) = TT^^^T^T^A^ (4.16) 



Then the traced bimetric equations of motion on maximally symmetric, proportional back- 
grounds become 

2(d-3)(ci-4) _ 
^ + (d_i)(d-2) - ^5 - ' (^-l^) 

2(d-3)(d-4) 2 .4^o^ 

^+ (d-l)(d-2)' ^2A-A/-0, (4.18) 

The above equations determine A and c in terms of the parameters of theory. For the 
purpose of PM analysis, we prefer that c is determined by a polynomial equation. For this, the 
solution for A can be easily obtained in the appropriate form on multiplying the /-equation 
by subtracting it from the g(-equation. This gives, 

A = — ^-^ . (4.19) 

C ^2 '2 

Then plugging this expression into the difference of (4.17) and (4.18) gives 
2((i-3)(d-4) 



{d-\)(d-2) 



■llilKg - c^AfY + (c^ - lc^){Af -Ag) = 0, (4.20) 



where I = 12/^2- This is the desired polynomial equation for c which is a generalization of 
the previously considered condition Af = Ag. In order to arrive at the PM candidate, we will 
demand that this equation leaves c undetermined. 

The effective Planck masses in terms of which the FP mass should be expressed now read, 



-"'3 V (^_l)(d-2)'2^V' V ('^-l)(^^-2) 

(4.21) 

These results are sufficient to investigate the most general PM theories in dimension d = 5 
and d = 6 where only the quadratic LL term can be added to the bimetric action. 
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4.3 The cubic Lanczos-Lovelock term 

The cubic LL invariant is given by 

+ SRR^^p^R^^'P'' + 24R'"'Rf'''R^p^„ + IQE^"^ Rf^R^p - URR^^R'"' , (4.22) 
which on maximahy symmetric backgrounds reduces to 

The trace of the equations of motion on proportional, maximally symmetric backgrounds 
then gives, 

2(d-3)(d-4) 4(d-3)(d-4)(d-5)(d-6) 3 _ 

^ + (d-l)(d-2) + (d-l)2(d-2)2 - - ' (4-24) 

2(d-3)(d-4) , . , 4(ci-3)(d-4)(d-5)(ci-6) , . 3 
^+ (d_i)(d-2)' ^2 A + (d_i)2(rf_2)2 ^ ^3 A -A/-0. (4.25) 

Subtracting off the cubic term will result in a quadratic equation for A and hence the solution 
will involve a square root. However, in order to apply our method for determining the PM 
candidate, we need to arrive at a polynomial equation in c. It turns out that this can still 
be achieved by decomposing the equations into a Grobner basis of equivalent polynomials 
instead. In the cases we have been able to reduce in this way it turns out that we then get a 
linear equation for A together with a consistency condition on the coefficients of the original 
polynomials which can be converted into a polynomial equation in c. Here, for simplicity, 
we present the results for the Grobner reduction after having imposed the symmetry relation 
In = a'^~'^"'ln which also covers the case of PM theories (although later the PM analysis is 
performed starting with the most general set of parameters). We then obtain a linear equation 
for A with the solution 

. ^ n3(l - g^)(l + q^A^A - nlq\l - q^)[\K + 2g^(Ag - Ay)] - n^n^j/^.kf 

{l-q2)[n^{l-q2){l+q^f-nlq^{l-q^-n2{kg-q'kf))] ' ^ - ^ 

where, we have used the following notion, 

d — 2ti 

AgA = Ag-q*Af, q = ac, nt = ^^Nk{d)ll . (4.27) 
A consistency condition on the coefficients of the original polynomials is, 

nl{Ag - q^Aff + nlq\l - q^) [q\Ag - A/ - naA^ ) - nsA^ - q^{Ag - A/ - 2n2AgAf)\ 
- n^q^{l - g2) [q\l + 3n2A/)(Ag - Af) + q\Ag - Af - naA^) - naA^ 

- q\l + 2,n2Ag){Ag -Af)- q\Ag -Af- 2n2AgAf)] = . (4.28) 

On multiplying by a factor c^^'^~'^\ one obtains a polynomial equation in c. The requirement 
that this leaves c undetermined constrains the /3„ and the In as given in Table 1. With these 
results for the cubic term we can investigate PM theories up to d = 8 in full generality. 
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4.4 Partial Masslessness with LL terms in d = 5 

Let us now consider the nonhnear PM candidate with higher derivative interactions in five 
dimensions where only the quadratic LL term is non-vanishing. This theory is parametrized 
by the six /3„, along with ^2 ^^"^ Planck masses tUg, tUf. In this case, (2.10) gives, 

5 

A9 = ^(/?o + 4c/3i + 6c2/32 + 4c3/33 + , (4.29) 
5 

A/ = ^(ac)-=^(c/3i + 4c2/32 + 6c=^/33 + 4c^/34 + c^/^s) • (4.30) 

Inserting these into (4.20) gives a polynomial with nine different powers in c. Demanding it 
to be independent of c by setting the coefficient of each power of c to zero therefore gives nine 
conditions on the parameters that can be explicitly solved. 

To simplify the presentation, here we first impose the extra requirement that the param- 
eters of the PM theory must correspond to a fixed point of the interchange symmetry (4.3), 
or else one ends up with two sets of nonlinear PM theories related by (4.3). Of course, we 
have checked that the unique nontrivial solution is also obtained without making this ansatz. 
Then this requirement gives, 

/3„ = a2"-'i/3,_^ , = a2-2";9 ^ (4.31) 

or, explicitly, /Js = a^/3o,/34 = a^/3i,/33 = 0/^2, /2 = a~^^2- Imposing these, it follows that the 
nine equations reduce to three independent ones, 

(a/3o - /3if - 3Ka/3i = , (4.32) 
2(a/3o - /3i)(a/3i - ^2) - ^nah = , (4.33) 
16(a/3i - 132? - l^Ko^h + 3Ka^/3o + 2,Ka^ (3i = , (4.34) 

where dimensionless parameter. The structure suggests the following definitions, 

/3o = — 3k6o , /3i = —Snabi , /32 = —3^0^62 . (4.35) 
Then the system reduces to the three equations, 

(60 - hf + 61 = 0, 2(60 - bi){bi - 62) + 62 = 0, 16(6i - 62)^ -bo-bi + 662 = 0, (4.36) 
which are solved by bo = ^,bi = —^,62 = ^, or, 

5 CK Q? Q? 

/34 = ^Ac, /35 = -^«, l{ = a-Hl. (4.37) 

These determine all the and I2 in terms of /f) "^g and mj. Inserting these into the solution 
for A in equation (4.19), and into the FP mass (4.14) gives the linear Higuchi bound, 

m\^ = ^{-\^ac- o?(?) = \^\. (4.38) 



- 19 - 



4.5 Partial masslessness with LL terms in d > 5 



In d = 6, the cubic LL term is topological so it is still enough to consider the quadratic term 
and the solution for the PM parameters can be obtained in a straightforward way. Up to a 
sign ambiguity in /33, all /3n are uniquely fixed in terms of /g- 

In d = 7, we also add the cubic Lovelock term and use the results of section 4.3. It turns 
out that there exists no nontrivial solution to the constraint equations on the /3„. Hence there 
is no PM theory in seven dimensions even with LL terms. 

In d = 8 it is still sufficient to consider only the quadratic and cubic terms. Remarkably, 
the resulting PM candidate theory has two free parameters. In particular, it is possible to set 
the cubic coupling to zero and obtain a PM theory with only the quadratic Lovelock term. 

For d = 9, 10 it has been checked that adding both £(2) arid /^(s), or adding £(4) alone, 
does not lead to nontrivial solutions for the parameters. The only possibility then is to add 
all terms up to quartic order, which we have not been able to investigate with our methods 
at hand. Hence we can only make definite statements about dimensions up to d = 8. 

The results for the PM parameter values are summarized in Table 1, where, to simplify 
the presentation, we define 

'^d - ig^d ' - ^igyi^d ■ ^^-^^^ 



Table 1. Parameters of the PM candidates 



dimension 


3 


4 


5 


6 


8 


'2 






1% 


l\ 


l\ 


li 






















1^ 
'3 


'3 












/3o 


/3o 


/3o 


12 '^S 


12 '^6 


"liso (20^8 ^ + 3Kg ^ j 


/3i 







12 ^^5 












2 


24 '^S 


12 '^e 


1600 '^8 


k 







24 '^S 







Pa 






12^*5 


12 '^e 


8000 l^^UKg yKg j 








5«5 (2) 
12 '^5 
















12 H 


1600 '^8 

























1600 (20«^8 + 3«^8 ) 
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Note that in all cases the parameters satisfy the symmetry property (4.31) and the linear 
Higuchi bound, 



d = 8: 



d = 6: 



d = 5: 



2 _ 



2 ^ 



TTlpp = - 




(4.40) 



5 Discussion 

Our results are summarized in the introduction. One of the issues left unanswered in the 
present analysis is the actual form of the gauge transformation in the nonlinear PM theory. 
Currently one can see the nonlinear form of this gauge transformation in the nonlinear theory 
in dS backgrounds only for constant gauge parameters. Non-constant gauge parameters will 
move the background away from dS. However, the analysis of cubic interactions [12, 14, 16] 
shows that, to cubic order, the theory is invariant under the full linear transformation, and 
not just its constant part employed in this paper. Another evidence for the existence of a 
nonlinear symmetry comes from the study of cosmological solutions in [36] based on a FLRW- 
type ansatz. It turns out that for the PM values of parameters (3.20) and the cosmological 
ansatz of [36] , the equations of motion leave a complete function of time undetermined which 
indicates a nonlinear gauge invariance. Furthermore, using a different approach, the authors 
in [15] consider partial masslessness in massive gravity in a certain decoupling limit and find 
a nonlinear symmetry that exists in that limit. This analysis is not directly comparable to 
construction here, but could be taken as another indication of the presence of the nonlinear 
symmetry. 
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A Mathematical details 

Here we provide the details for deriving the equations of motion from the bimetric action 
(2.1), computing the cosmological constants on proportional backgrounds and finding the 
Fierz-Pauli mass of the linear fluctuations. 



A.l Bimetric equations of motion 



The bimetric interaction potential in (2.1) is given in terms of en{S) where S = \/g\f- For 
a d X d matrix S, one has 

d 

det(l + A5) = ^A"e„(5). (A.l) 



n=0 



e„(5) are the elementary symmetric polynomials of the eigenvalues of of S and can be itera- 
tively constructed starting with 60(5*) = 1 and using the Newton's identities, 



n(S) = -if;(-l)'=TY(5^-)e„„fe(5) 



(A.2) 



k=l 



In d dimensions, the iteration ends with ed{S) = det S and then e„(S') = for n > d. 
Obviously, en(AS') = A"'e„(S'). To obtain the equations of motion, one needs the variations. 



5en{S) = - ^(-l)"Tr(S'"-i55) e„_^(5) . 

m=l 

These follow from (A.2). Hence, 

n 

5[dei{S-^) en{S)] = - det{S-^) (-l)™Tr(5"^-i^5) e„_„(5) . 



m=0 



If S is the square-root of a matrix E, = E, then 

1, 



Tr(5'"-V5) = -Tv{^/E 5E) 



(A.3) 



(A.4) 



(A.5) 



Note that for matrices, this relation holds only under the trace. If E = g then on varying 
g~^, changing the summation variable from m to r = n — m, and cancelling y^|7|, one gets. 



where, the matrices Y^^^[S) are given by. 



(A.6) 



(A.7) 



fc=0 



Using (A.6), for V given by (2.3) and satisfying the property (2.5), it is now straightforward 
to compute, 

d-l 



V3 



-W\9\V 



(A. 



n=0 



I/I Sf 



-2V\T\v] = Ux j;(-i)"/3,_„y(t),(v7^) • (A.9) 



n=0 



This leads to the bimetric equations of motion (2.6). Note that the two equations are obtain- 
able from each other by the interchanges g^i, -H- f^^, f^n ^ f^d-n, f^g -H- mj. 
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A. 2 Proportional backgrounds and the cosmological constants 

Now consider the equations of motion (2.6) for the background ansatz f^^ = c^g^,^, implying 
S = cl. It is obvious that the potential terms (A. 8), (A. 9) become cosmological contributions, 

-^V^, = g,,Ag, -^vf, = g,,Af. (A.IO) 

Ag and Af are expressed in terms of y(„)(c]l). To compute this, consider (A.l) for S = cl. 
Since det(l + XS) = (1 + Xc^ and efc(cl) = 6^(1), it follows that 

Exactly this sum appears in Y(„)(cl) (A. 7) which becomes, 

y(„)(cl) = (-l)"c"(^^- ^) , F(„)(c-il) = (-l)«c--(^^- ^) . (A.12) 
On substituting in V^^^ one reads off the cosmological constants Ag and Af from (A.IO) as, 

rrip \ n J rUf ^ \n-lj 

y n=0 / n=l 

Therefore, for proportional backgrounds, the bimetric equations reduce to two copies of Ein- 
stein's equations (2.9). The consistency of these background equations with each other then 
requires A^ = Af, which generically determines c. 

A. 3 Fluctuations and the Fierz-Pauli mass 

Let us now consider linear perturbations around the proportional backgrounds (2.8), 

9tiu = (jfiu + Sg^u , f^lu = Uu + 5f^y , (A.14) 

with / = (?(}. The corresponding linearized equations obtained from (2.6) are, 

6 {R,,{g) - y^^uRig)) + ^ 6V^, = , 6 [R^,{f) - \U.R{f)) + SV^, = . 

(A.15) 

To compute these explicitly, we need, 

^^^A-S = ^A,5ff,, Y.^-lr^r^5Y^r^)u{S)\-s (A.16) 

n=l 

,d-2 d~l 



Till 



n=l 



The first term, obviously, is the contribution of the cosmological constants to the linearized 
fluctuation equations. The second term contains the FP mass and involve 5g and 6f only 
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in the combination SS'^^ = ^g^^{Sf — (?6g)\y. Hence the massive mode must contain this 
combination of the fields. Since 5Yq = 0, the mass term is independent of both /3o and /3rf. 
From (A.7), 



r=0 

The first term in (51^(^)1, proportional to 5S, involves, 

n 

nJ^(-ire.(l)=n(-irC^^ ^) =(d-l)(-ir 

r=0 

n 

Y,{-lYrer{l) = d{-l) 



(A.18) 



r=l 



n 
d-2 



d-2 
n — 1 



n — 1 



(A.19) 
(A.20) 



To compute the second term in (51^(^)1, proportional to Tr(5S'), note that (A. 3) and (A. 11^ 
give, 



5er\ 



: -c'-' Tr(<55) j2 (-1)" ^r-mil) = c'~' Tr(<55) _ J) • (A.21) 
with (5eo = 0. Then, the second term in (A.18) contains, using the second equation in (A. 11), 



Putting these together in (A.18) gives, for n > 1, 

'd-2 



5y(„)(5)|^ = (-i)V 



n n— 1 



n — 1 



[Tr(55)l - 6S] . 



(A.23) 



From this one can easily obtain 5Y(„)(S' ^)|^ by replacing c — 1/c and noting that 5{S ^] 
—c~^5S. Finally, putting all this together one gets. 



6V^ - 



5V^ I - 



m„ 



,d-2 



m. 



Afc--'6f^x - c^-^Ng^x {"^^{53)5^ - S3 



with. 



N 



.71=1 ^ ^ 



(A.24) 
(A.25) 

(A.26) 



The Einstein tensor in (2.6) is linearized in the standard way: one has, 



S ( R^u + -j^gfivR ] = SjluSgpa + 'R-tiu ■ 



(A.27) 



- 24 - 



The structure of £{luSgpa is given by (1.2) but with the background metric g^j^, and, 

= -\ {Mg^v - %uR"'^9pa) = {dbg^, - g^ug'^Sgp^) (A.28) 

In the last step we have used the background equation (2.9). Finahy, using (A. 27) and the 
corresponding equation for 5/, along with (A. 24) and (A. 25) in (A. 15), we arrive at the 
linearized equations (2.13) and (2.14) for the fluctuations 5g and 5f. 

B The nonlinear G-M'^ action 



We sketch a brief derivation of the nonlinear bimetric action (2.1) written in terms of the 
fields 

G^, = (l + {acf-^) ^'"'^/^'^'^ (^g^^ + {acf-^c-^f^,) , (B.l) 

and 

^ / J ^\ -2/(d-2) 

Mg. = (l + {acf-^) {G^,SP, -cG^,) , (B.2) 

where we recall that S = sj g~^f. These definitions are consistent with all of our consid- 
erations in the main text. From these definitions we straightforwardly obtain the inverted 
relations 

g,. = G,p{^-y, , V = G,aS\S\{^~y, , (B.3) 
where, for notational purposes, we have defined 

= (l + (ac)^-2\ (4-'i)/(<^-2) / ^ ^ (^acr-2,-2sP^s%) . (b.4) 



From (B.2) we can also directly obtain 5 in terms of G^^ and 



SP^ = c6i:+{l + {ac)^-') GP>^m;^^. (B.5) 

It is now easy to see that the volume densities can be rewritten through 

V5 = \/Gdet($)-^/2^ v7 = \/Gdet(S)det($)-i/2_ (gg) 

In order to rewrite the interaction potential we note a general relation obeyed by the sym- 
metric polynomials, 

e„(l + X) = ^f _^ e,(X). (B.7) 
fc=0 ' 

If we consider the matrix parameterization 5 = c (1 -|- X) (c.f. (B.5)), we may use this relation 
to rewrite the interaction potential, noting that 

d d 

V[S- fin) = /3ne„(5) = c"«„e„(X) = y(X; c^a^ , (B.8) 

n=0 n=0 
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where the different set of coefficients are related by 

d / , \ d 



c an = 

k=n ^ k=n 



This is straightforward to implement in the interaction potential, which is then given by 

V{S;(3n) = V{G-^M^;an), = (l + {acf-^j a„ . (B.IO) 

For the kinetic terms we choose to write these in terms of covariant derivatives with respect 
to G^^. Using the Riemann curvature definition [V^,V,y]u;p = Rf^^p^Ua, together with the 
general relation between two covariant derivations on a manifold, 

V^^. = V^:^. - C^/ujp , Cp/ = i/'^ (2V«(^5,), - V^/^.) , (B.ll) 

we obtain the relation (with an obvious similar relation for Ri_iuif) obtained by replacing 
g — 7- / in these expressions) 

R,u{g) = RAG) - 2V£C ]/ + 2C [/C ]/ . (B.12) 

Tracing this with g^^'^ {f^^) and writing it out in full we find that (modulo total derivatives 
and neglecting to write out an overall factor of ■\/\g\ (\/[7I) both sides) we get Ricci 
curvature relations of the form 

R{g) =g^'RAG) - i^pAV^^^^V^^"" + \gp.^ xg"""^ aO"'' 

- \gp.gxpg''''yag'^yag^^ + iff./3<7pAff"'^V,9"^V./^ ^ g'^^'RAG) + , (B.13) 



with a similar relation for R{f) and definition of II-^, obtained by replacing g ^ f in this. 

i^u and M^^ 



Collecting our results we can now write the nonlinear action in terms of Gnu and Af^ as 



Sgm = 2 J d^:,^\det{<^>r'/^ [{^G-'rR^G) + a^~^ det{S){<^G-' S~^r R^G) 

+ + det(5) - 2^ V{G~^M^; On)] , (B.14) 



where the relations (B.3), (B.4), (B.5) and (B.13) can be used to get the explicit form in 
terms of only and . Expanding this action to second order gives the quadratic action 
(3.25) by construction. 
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